Effective potential at finite temperature in a constant hypermagnetic field: Ring 

diagrams in the Standard Model 
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I. INTRODUCTION 

The problem of baryogenesis is still one of the out- 
standing open questions in cosmology, despite the large 
amount of work devoted to find a viable explanation. The 
conditions for developing a baryon asymmetry in an ini- 
tially symmetric universe were established by Sakharov 
in 1967 [H and the search for a scenario to encompass 
them continues. These three well-known conditions are: 

(1) existence of interactions that violate baryon number; 

(2) C and CP violation and (3) departure from thermal 
equilibrium. 

The Standard Model (SM) of electroweak interactions 
meets all these requirements, provided the electroweak 
phase transition (EWPT) be first order, since, at that 
stage of the universe evolution, this is the only possible 
source of departure from thermal equilibrium. Nonethe- 
less, it is well known that neither the amount of CP vio- 
lation within the minimal SM [2| , nor the strength of the 
EWPT are enough to generate a sizable baryon num- 
ber [|. 

During the EWPT, the Higgs field vacuum expectation 
value changes from zero (false vacuum) in the symmetric 
phase, to a finite value (w) (true vacuum) in the broken 
phase. This evolution is determined by the finite tem- 
perature effective potential of the theory, that develops 
a barrier between the two minima if the phase transition 
is first order. The temperature at which the two min- 
ima are degenerate is called the critical temperature 
and this instant is considered the beginning of the phase 
transition. From there on, the transition is accomplished 
through the nucleation, expansion and percolation of true 
vacuum bubbles in the background of false vacuum, lead- 
ing to a departure from equilibrium conditions. 

The existence of baryon number violation is realized 
in the SM by means of its vacuum structure through 
sphaleron mediated processes. The sphaleron is a 
static and unstable solution of the field equations of 
non-Abelian gauge theories, corresponding to the top of 
the energy barrier between topologically distinct vacua, 
where the minima correspond to configurations with zero 
gauge field energy but different baryon number. Transi- 



tions are associated to baryon number ns — violation 
and can either induce or wash out a baryon asymmetry. 

The above property of sphalerons makes the preserva- 
tion of a given baryon asymmetry one of the most difficult 
conditions to meet during the baryogenesis process in the 
majority of the proposed scenarios. This requires that 
the baryon violating transition between different topo- 
logical vacua is suppressed in the broken phase, when the 
universe returns to thermal equilibrium. In other words, 
the sphaleron transitions must be slow compared to the 
expansion rate of the universe and this in turn translates 
into the condition (v)/Tc > 1.0 - 1.5 Q. Although the 
above is an estimate emerging from approximate calcula- 
tions, it is nowadays widely accepted and has proven to 
be a rather difficult condition to meet. 

Although there have been several attempts to link the 
baryogenesis process to the EWPT 0, Q in general these 
all share the characteristic that the Sakharov conditions 
are only partially met. Here we want to further ex- 
plore the possibility to embed the baryogenesis process 
in the EWPT scenario, including the effect of an extra 
ingredient: the possible presence of primordial magnetic 
fields. Before the EWPT, magnetic fields couple to mat- 
ter through the particle's hypercharge and thus properly 
receive the name of hypermagnetic fields. 

Magnetic fields seem to pervade the entire universe 
and their generation may be either primordial or asso- 
ciated to the process of structure formation. They have 
been observed in galaxies, clusters, intracluster medium 
and high redshift objects In order to distinguish 
between primordial and protogalactic fields, it is use- 
ful to search for their imprint on the cosmic microwave 
background radiation (CMBR). A homogeneous mag- 
netic field would give rise to a dipole anisotropy in the 
background radiation, on this basis. Cosmic Background 
Explorer (COBE) results give an upper bound on the 
present equivalent field strength of Bq < 10"^ G 8]. An 
upper limit to the field strength of Bq ^ 4.7 x 10^^ G 
at the present scale of 1 Mpc is obtained by an anal- 
ysis that includes small scale CMBR anisotropics from 
Wilkinson Microwave Anisotropy Probe (WMAP), Cos- 
mic Background Imager (CBI) and Arc Minute Cosmol- 
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ogy Bolometer Array Receiver (ACBAR) On the 

other hand, tangled random fields on small scales could 
reach up to ^ 10~^ G Nucleosynthesis also imposes 
limits on primordial magnetic fields since they have an 
influence on both the universe expansion rate and the 
electron quantum statistics. The observed helium abun- 
dance implies B ^ 10^^ G at scales greater than 10 cm 
at the end of nucleosynthesis [ll|. 

Although at present there is no conclusive evidence 
about the origin of magnetic fields, their existence prior 
to the EWPT cannot certainly be ruled out making it 
important to investigate their effect on the baryogenesis 
process [l^. In fact, it has been shown that these fields 
provide mechanisms to affect all the Sakharov conditions: 
In the presence of a magnetic field, the phase transition 
becomes stronger first order in analogy to the case of a 
type I superconductor, where the Meissner effect brings 
the phase transition from second to first order [l3l.[T3.[l5|: 
it has also been shown that extra CP violation is obtained 
from the segregation of axial charge during the reflec- 
tion and transmission of fermions through the vacuum 
bubbles due to the chiral nature of their coupling to the 
hypermagnetic field [l6| ; finally, regarding the sphaleron 
transitions the presence of a magnetic field works against 
the preservation of a baryon asymmetry due to the cou- 
pling between the sphaleron dipole moment and the mag- 
netic field that lowers the energy barrier between topo- 
logically distinct minima ■ 

The effect of magnetic fields on the EWPT has been 
analytically studied both classically [l^ and to one- loop 
order dl| , as well as by means of lattice simulations [l5| . 
These calculations all agree that the strength of the phase 
transition is enhanced by the presence of hypermagnetic 
fields, although the ratio (v) /Tc does not reach the de- 
sired value, for a large Higgs boson mass. On the other 
hand, other analytical approaches where the SM finite 
temperature effective potential is studied for the case of 
strong magnetic fields [l^ [l^ , reach the conclusion that 
these fields inhibit the first order phase transition and 
attribute the result to the contribution of light fermion 
masses which are generally neglected in other computa- 
tions. 

In this work we concentrate on studying the relation 
between the presence of a large scale magnetic field and 
the dynamics of the EWPT by computing the SM finite 
temperature effective potential in a constant hypermag- 
netic field up to the contribution of ring diagrams, that 
have been shown to be crucial for the des crip tion of the 
long wavelength properties of the theory [2(l|. We carry 
out a systematic calculation of each SM sector showing 
that the major contribution producing an enhancement 
of the EWPT comes from the Higgs and gauge boson sec- 
tors and that fermions do not act against this behavior. 
Working in the limit where the magnetic field is weak, 
we find an enhanced value of the ratio (v) /Tc. 

The paper is organized as follows: In Sec. |TT] we lay 
down the formalism to include weak magnetic fields in 
the computation of hypercharged particle propagators. 



In Sec. mil we write down the SM using the degrees of 
freedom in the symmetric phase. In Sec. IIVI we work 
with these degrees of freedom to compute particle self 
energies that are used in Sec. |V] to compute the SM ef- 
fective potential up to the contributions of ring diagrams. 
In Sec. IVII we study this effective potential as a function 
of the Higgs vacuum expectation value and show that the 
order of the EWPT becomes stronger first order in the 
presence of the hypermagnetic field. Finally, we conclude 
and discuss our results in Sec. IVIII and leave for the ap- 
pendix the computation of some intermediate results of 
the analysis. 



II. CHARGED PARTICLE PROPAGATORS IN 
THE PRESENCE OF A HYPERMAGNETIC 
FIELD 

We work with the degrees of freedom of the SM in the 
symmetric phase, where the external (hyper)magnetic 
field belongs to the U{1)y group. To include the ef- 
fect of the external field, we use Schwinger's proper time 
method '^l*]. In the symmetric phase, we have only two 
kinds of hypercharged particles that couple to the exter- 
nal field namely, scalars and fermions, whose propagators 
are 



Dh{x,x') — 0(x,a;') 



Sh{x, x') — 4>{x, x') 



Dnik), (1) 



(2^ 



«fe.(a,-x')5'^(fc)^ (2) 



respectively. The phase factor (j){x, x'), that breaks trans- 
lation invariance, is given by 



{x,x') 



(3) 



where the vector potential -B™' — (0; ?/j 0) gives rise 
to a constant hypermagnetic field of strength H along the 
z axis and = d^B'i^^ — d^B°^^ is the external field 
strength tensor. 

The momentum dependent functions Dnik) and 
Snik) are given by 



tDHik) = / 
Jo 



cos Y_ff S 
2 



X exp <{ is{kf, ~ k\ - + ie) \ ,(4) 



iSnik) 



ds 



cos yHs 



exp |js(fc^ — k^ 
{ruf - |(f||)e' 



tan yHs 



m + ie] 



cos yTJs 



(5) 
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where y is the particle's hypercharge and we use the no- 
tation /cj^ = /cq — /e| and fc^ = fcj + fc| . 

Since the gauge bosons do not couple to the external 
field their propagator is given by 



9' 



(1-0 



k^-^m'i^ 



fc2- 



(6) 



where a,b = 1,2,3,4 and the first three values correspond 
to the SU{2)l fields and the fourth to the U{1)y field. 
Notice that the matrix m^. is not diagonal in the basis 
of the weak-interacting fields in the symmetric phase. 

It has been show that, by deforming the contour of 
integration, Eqs. ^ and ([S]) can be written as [13, 



iDnik) = 2i 



(-iyL,(^)exp{-^} 



{21 + 1)yH - m2 



(7) 



^ k?, - 21yH - tri^f + ie 



1=0 "-11 



where di{a) = {-l)-^e-~°'L^\2a), d'^ = ddn/da, 



is obeyed, were we consider m as a generic mass of the 
problem at the electroweak scale. 

We can thus perform a weak field expansion in Eqs. ([7]) 
and ([5]) which allows to carry out the summation over 
Landau levels to write the scalar and fermion propagators 
as power series in 'yH, that up to order (yH)^ read as (2^ . 

m 



D{k)H 



1 



fc2 



(Yg)2 2{^Hfkl 
(fc2 — m2)2 (fc2 — m?Y 



(11) 



and 



S{k)H 



ji + mf 75^^(fc|| -Hm/)(Yi/) 



2\3 



f-mj {k^-m}) 



(12) 



respectively. 

There is an analogous result for gauge bosons, but since 
in the symmetric phase these do not couple to the exter- 
nal hypermagnetic field, we do not need to account for 
them. 



mi — /cm 



D 



IbiPimf + fill), 



(9) 



Li , Lf^ are Laguerre and Associated Laguerre polynomi- 
als, respectively, and u'^, b'^ are four- vectors describing 
the plasma rest frame and the direction of the hyper- 
magnetic field, respectively. 

In order to set the appropriate hierarchy of energy 
scales we resort to qualitative cosmological bounds on 
the possible strength of primordial magnetic fields dur- 
ing the EWPT. CMBR sets stringent bounds on large 
scale primordial fields but no so much stringent when the 
fields are tangled. This dependence on the scale makes it 
difficult to extrapolate these bounds down to the EWPT 
epoch. We use instead the requirement that the magnetic 
energy density Pmag ~ should be smaller than the 
overall radiation energy density prad ~ T'^ at nucleosyn- 
thesis, in order to preserve the estimated abundances of 
light elements. With this, one obtains the simple bound 

On the other hand stability conditions against the 
formation of VF-condensate indicate [25] that the field 
strength is also weak compared to the square of the 



W mass, 



Notice however that when thermal cor- 



rections are taken into account, this bound could be 
avoided 

We work explicitly with the assumption that the hier- 
archy of scales 



YiJ < TO^ ^ 



(10) 



III. STANDARD MODEL 

In order to consider all the contributions to the SM 
effective potential, we write the Lagrangian for each sec- 
tor. 

The Lagrangian for the Higgs sector is 
Ch = {D^<S>)\D''<i>) + c2($t$) _ A($^$)2, (13) 

where = d,, + ig^A^^^ + i^B'^, r° are the Pauh 
matrices, B'^ = + and B^ are the SU{2)l 
and U{1)y gauge bosons, respectively. To allow for spon- 
taneous symmetry breaking the mass parameter c must 
satisfy > 0. 

The Higgs field is a complex doublet with y = +1 



(14) 



where ipi are real scalar fields. We take (j>i as the physical 
Higgs field that develops a vacuum expectation value v. 
Extremizing the tree level potential, the parameter c is 
related to the classical minimum Vdass by 



2 \ 2 

c = Aw, 



class • 



(15) 



For the Higgs field hypercharge conjugate doublet we 
use $ — icr2$*. 

The kinetic energy from the SU{2)l and U{1)y gauge 
bosons is 



C 



1. 



gb 



1 pffiiy pf 

4 



(16) 



where 

K'^ = d^.Bl~d.B'^. (17) 
The Lagrangian for the fermion sector is 

(^z^-|^Y^'-|r.A^*i, (18) 

where *l = i(l - 75)^' is an SU{2) doublet of left- 
handed fermions and 'i'u — |(1+75)'J' is an SU{2) singlet 
of right-handed fermions. 

We work in the limit that all fermion masses except 
the top quark mass are negligible, so that the main con- 
tribution to the Yukawa sector is 

^Yukawa = fW^tp. + h.C. (19) 

Finally, in the i?^ gauge, the gauge fixing Lagrangian 

is 

Y^{d^B,-\Uv<t>2)\ (20) 

where « = 2, 3, 4 and ^ is the gauge parameter. We choose 
to work in the Landau gauge (^ = 0) in which the ghost 
fields do not acquire mass and hence do not contribute to 
the D-dependent part of the one-loop effective potential. 
Note that the effective potential is in principle a gauge 
dependent object [HI, however, physical quantities ob- 
tained from it are gauge independent [27j . 
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FIG. 1: Self-energy Feynman diagrams for the Higgs bosons 
that contain loop particles affected by the hypermagnetic 
field. These particles are represented by thick lines. "1> and 

represent Higgs and Fermion fields whereas AJJ and 
represent the U(V)y and SU(2)l gauge fields, respectively. 




FIG. 2: Self-energy Feynman diagrams for the Higgs bosons 
that contain loop particles not affected by the hypermagnetic 
field. These particles are represented by thin lines, c repre- 
sents the ghost fields. Working in the Landau gauge ^ = 0, 
the second diagram vanishes. 

A. Higgs boson 



IV. SELF-ENERGIES 

In this section we compute the SM self-energies that 
are in turn used for the computation of the ring diagrams 
in the effective potential. 

It is well known that in the absence of an external mag- 
netic field, the SM thermal self-energies are gauge inde- 
pendent when considering only the leading contributions 
in temperature [28]. However, as we will show, when con- 
sidering the effects of and external magnetic field, these 
self-energies turn out to be gauge dependent. 

In what follows, we work in the imaginary-time for- 
malism of thermal field theory. First, we note that the 
integration over four-momenta is carried out in Euclidean 
space with fco = ik^, this means that 

Next, we recall that boson energies take discrete values, 
namely A:4 = a;„ = 2mrT with n an integer, and thus 




Figure [T] shows the diagrams that contribute to the 
Higgs boson self-energies affected by the hypermagnetic 
field. Let us explicitly compute the momentum indepen- 
dent diagram shown in fig. [Ij/) for a single scalar field. 
In the weak field limit, its expression is 

W "4^7 (27r)3 

X DH{uJn,k;m^ ^ (23) 

where Dh is given by Eq. pT|) . Notice that Eq. ([23|) is 
computed self-consistently, with the approximation that 
on the right hand-side, H^^^^* Hi, where Hi is given 

by m 

ni = ^{l9' + \g" + 2X + fY (24) 

and represents the leading temperature contribution to 
the scalar self-energy. The need to compute H^*^^'* self- 
consistently is linked to the fact that, on the one hand, in 
the SM, scalar masses can vanish as a function of v and, 
on the other, the presence of the magnetic field origi- 
nates terms inversely proportional to these masses [see 
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Eq. p2|) ]. Thus, for soft momentum, where the contribu- 
tion of the ring diagrams is relevant, a naive perturbative 
expansion is not sufficient. The well known correction of 
the infrared behavior is given by the plasma screening 
properties and in the case of Eq. (|23p we approximate 
such correction as consisting only of the leading temper- 
ature contribution to fl^^f^'* which, upon resummation, 
take care of the most severe infrared divergences 26] . 
Using the Euclidean version of Eq. (fTTj) , we have 



n 



Higgs 
(!) 



A 



1 



(27r)3 (^2 + k2 m2) 

(YiJ)2 



2 
± 



2(Y7J)2k2 



(tj2 -(- -I- rh?Y 



(25) 



where we use the short hand notation — rri^ -(- Hi . 

The integrand in Eq. (j25p contains terms whose general 
form is 



/a/3(k,q) = 



1 



(26) 



We make use of the Feynman parametrization to write 

Ial3 as 



where 



Via + (3) dx x'^-^il- x)!^'^ 



r{a)T{(3) Jo [uol + \<i'\x) + m'2(x)]"+^ ' 



(27) 



k'(a;) = k-(l-a;)q 



to' {x) 



a;(l -x)q^ 



(28) 



and r is the Gamma function. Notice that this 
parametrization is allowed since we work in the 
imaginary-time formalism (29j . 

To carry out the sum over Matsubara frequencies to- 
gether with the integration in Eq. (|25p . we perform 
an asymptotic expansion in the high temperature limit. 
This is done by means of a Mellin transform as described 
in detail in Ref. [s^. The explicit result, generalized to 
also include the fermion case is 



TV / ^'^k2a,,2*, f,. (2r)(2.r)'^+2.+2t-2(.+,)r(|+a) 



(27r)° 



(47r)'^/2r(a)r(/?) 



r(f) 



X4^C(2(j- + « + /?-t-^-a),Z) 



J'- 



X r(j + a + /3-^-a) fdxx"-^ (i _ ^)^-i ("l^ 



2j 



(29) 



r 



where C is the modified Riemann Zcta function, fj, is the Eq. ((29|) calls for the use of dimensional regularization, 

energy scale of dimensional regularization and d = 3 — which is an appropriate method to use for non-Abelian 

2e. For fermions the sum runs over all integers n and gauge theories. 
Z = 1/2, while for bosons the n = term is excluded 

and Z ~ 0. It is important to stress that the method For the terms involving the n = Matsubara frequency 

advocated in Ref. [30] to perform an expression such as for bosons, we use the result 



T I ^k-/ .(k q) - r r(i-ha)r(a + /3-f) _ 

^ J {2n^^f ^-^^^^"^^ - [A^yn r(f) r(a)r(/3) X^^^ \m'{x) 



2a+2f3-d+2a 



(30) 



In terms of the functions Iaf3 defined in Eq. (|26p we 
can write Eq. (p5)) as 



/4o(k,0) 



(31) 



n 



Higgs ^ ^rp\^ 



{I) 



A ^ J (27r)3 



/io(k,0)+/3o(k,0) 



Using Eq. for the terms with n 7^ and Eq. 
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for the term with n = 0, we get 



TjHiggs _ ^rp2 

(I) 2 



2 

E 

1=1 



3 TO, (YiJ)2 

^ IGTrTm? 



(32) 



where we have included the contribution from all scalar 
fields with rfif = mf + IIi, to^ standing for the scalar 
boson masses, given by 



4° R ^ 




Too 



3Aw' 

™2 _ 



(33) 



In a similar fashion the contribution from the diagram in 
fig.[T]^//) in the infrared limit, namely go = 0, q ^ 0, is 



n 



Higgs 
{II) 



(0) 



1 



14C(3) {YHf 



{2n) 



rp4 



(34) 



where hereafter we use the notation J-{0) = J-{qo = 
0, q — > 0) to represent the infrared limit of any function 
J'. 

We point out that in Eqs. ([5^ and we have 

kept terms representing the leading contribution of each 
kind arising in the calculation, namely, terms of order 
{^Hf/T^, m,/T and {y:Hf /Tm^^. For the hierarchy of 
scales considered, the first kind of terms can be safely 
neglected. Recall that terms of order fhi/T are usually 
neglected in a high temperature expansion. However, 
since we are here interested in keeping the leading con- 
tribution in the magnetic field strength, we are forced 
to keep this kind of terms which, for a large top quark 
mass, namely a large /, are of the same order as terms 
{'yHY /Tfh^. Notice that for large values of the coupling 
constants, a perturbative calculation is not entirely justi- 
fied. Nevertheless, here we consider our calculation as an 
analytical tool to explore this non-perturbative domain 
and regard terms of order ifii/T as an estimate of the 
theoretical uncertainty of our results. 

On the other hand, the diagram in fig. [Ij///) is pro- 
portional to the parameter ^ and thus vanishes for our 
gauge parameter choice. Accounting also for the dia- 
grams that are not affected by the hypermagnetic field, 
and that are depicted in fig.IH the Higgs field self-energy, 
in the infrared limit is given by 



n^*99^(0) = 



' ^3(TOf+Hi)l/2 



T 



8T(TO2-f ni)3/2 



(35) 



^AAAAAAA/\?^y^7vWWWVV" 
(III) 



FIG. 3: Self-energy Feynman diagrams for the gauge bosons 
that contain loop particles affected by the hypermagnetic 
field. These particles are represented by thick lines. "I> and 
'5 represent Higgs and Fermion fields whereas A^^ and 
represent the U{1)y and SU{2)l gauge fields, respectively. 



to the gauge boson momentum g'^, namely u^, q^^ and 
bf^. This means that in general, these self-energies can 
be written as linear combinations of nine independent 
structures [3l| . Since we are interested in considering 
the infrared limit, go = 0, q — > 0, only u'^ and 6^ re- 
main. Notice that the correct symmetry property for the 
self-energy is H^^(9) ~ ^abi~l) I^ll- However, in the 
infrared limit, this condition means that the self-energy 
must be symmetric under the exchange of the Lorentz 
indices and therefore we can write. 



M LIV 



(36) 



where 



(37) 



and the transversality condition g^^H^^^ = is trivially 
satisfied in the infrared limit. 

Figure [3] shows the gauge boson self-energy diagrams 
that are affected by the external hypermagnetic field. 
These include diagrams involving scalars as well as 
fermions in the loop. Let us explicitly compute the dia- 
gram shown in fig. ^I) for the field. Its expression 
is 



_(£k_ 



DH{k)DH{k-q). 



(38) 



B. Gauge bosons 

To express the gauge boson self-energies, in the pres- 
ence of the external field, we have three independent vec- 
tors to our disposal to form tensor structures transverse 



Notice that since the net hypercharge fiowing in the loop 
is zero, the phase factor in Eq. ^ vanishes. 

Let us compute 11°°^^ for a single scalar field with mass 
TO. Working in the rest frame of the medium, n(/)B ~ 



H' 



{I}B 



ttM 

^hi)- 



From Eq. 



and at finite temperature 
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FIG. 4: Self-energy Feynman diagrams for the gauge bosons 
that contain loop particles not affected by the hypermagnetic 
field. These particles are represented by thin lines, c repre- 
sents the ghost fields. 



this component is given by 



n 



00 

(I)B 



(Pk 



^ J (27r)3"'" 

X Djy(w„,k;TO^) 
X I?//(tj„,k- q;m2). 



(39) 



Since the integral in Eq. pQ]) does not give rise to terms 
inversely proportional to m, we have omitted the replace- 



ment m? 



Hi. 



Using the Euclidean version of the scalar propagator 
Dh obtained from Eq. (jlip , and in terms of the functions 
/a/3 defined in Eq. ^7\ . we can write Eq. ([55)1 as 



J 



^"(1)3(0) = ~9"tY^ J ^c^?.{/ii(k,0) + (Yi/)2[/i3(k,0)+/3i(k,0) + 2ki/i4(k,0) + 2kiJ4i(k,0)]}. (40) 



r 



Using Eq. ([^ into Eq. (|^ and keeping only the lead- 
ing term as discussed in Sec. IIV Al we get 



n 



?zV(0) = 



(41) 



where the contribution from the four real scalar fields has 
been accounted for. 

In a similar fashion, the explicit expressions for n9jj-,g 



and ^^jjj^s are computed to yield 



^tl 1)3(0) 



n' 



{III)B 



(0) 




where in the first of the Eqs. (jH]) we have performed the 
sum over all hypercharged fermions. We stress that the 
origin of the terms ~ 1 /fhf is the topology of diagrams 
such as the one in fig. ^III) or fig. [ij/), involving a 
tadpole of hypercharged scalars in the presence of the 
external field. In the computation of these diagrams, we 
require to consider the replacement —^ m? + Hi (see 
the discussion in Sec. IIV A|) . 

Adding up the above three contributions, we get 



6 ^ 



1 - 



4 r 

E 

i=0 



3 fhi 



1 {^Hf 



7047r Tmj 



(43) 



Performing a similar exercise for the case of the 
fields, for which we also include the contribution from 
the diagrams that are not affected by the hypermagnetic 
field, depicted in fig. SI we can write the result for the 
four gauge bosons as 



ttOO 



,11 , 




3 (mf-fHi 



,1/2 



887r 



T 



7047rT(m2-hni)3/2 



Sab 



(44) 



where g = g for a=b=l,2,3 and g = g' for a=b=4. 

As is sketched in the appendix, the other non zero 
components of the gauge boson self-energy are negligible 



T-rll _ xt22 tt33 

^^ab — ^^abi ^^afc 



n 



n 



(45) 



which means that 11°° 



n 



ab ' 



Notice that in Eq. (M)) we have not included terms 
of order 0{Mab/T), where Mab is the gauge boson mass 
matrix. These masses are proportional to gv whereas for 
large /, •fhi are proportional to fv. Thus in this limit, 
terms of order 0{Mab/T) are smaller than terms of order 
0(jhi/T). These terms can in principle be calculated by 
diagonalization of Mab or else by explicitly considering 
the calculation in the basis of fields in the broken sym- 
metry phase. We will present such calculation in a forth- 
coming work along with the gauge parameter dependence 
of the self-energies and effective potential. 
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Although in principle the fermion self-energies are also 
affected by the hypermagnetic field, as in the case of zero 
external field, their contribution to the ring diagrams is 
subdominant in the infrared and do not need to be taken 
into account. 



where in general each one of these contributions is given 
by 



EFFECTIVE POTENTIAL 
A. One-loop 



y'''(-)-^T. J (03Tr(ln[i?K,k)-]), (48) 



In the standard model the tree level potential is 



(46) 



To one loop, the effective potential (EP) receives contri- 
butions from each sector, namely 



(47) 



with D stands for cither the scalar, fermion or gauge 
boson propagator, and the trace is taken over all internal 
indices. 



In the weak field limit, the contribution from the Higgs 
sector is given by 



V, 



(1) 



Higgs 



E I E / (03 HD,' (u.„, k; - + n,)] - E I E / (03 { H^l + + + n, 



i—1 n 
\2 



i—1 71 



2{kl) 



} 



r 



(49) 



The first term in Eq. (|^^ with Hi = represents the low- 
est order contribution to the effective potential at finite 
temperature and zero external magnetic field, usually re- 
ferred to as the boson ideal gas contribution [28| . In order 



to keep track of the lowest order corrections in A, we set 
Hi = in Eq. . Thus for the hierarchy of scales con- 
sidered here and dropping out the zero-point energy, this 
contribution is given by [26| 



J 



E|E/7gjM-^+k^+-?)-E(- 

1 n V / ,-—1 V 



90 



mjT^ mfT 



24 



In (- 



12n 327r2 V47rT 



+ Oimt) 



(50) 



r 



Notice that there are potentially dangerous terms mf in 
Eq (|50p that can become imaginary for negative values 
of rrii . However as we will show, these terms cancel when 
including the Higgs contribution from the ring diagrams. 
The second, if-dependent term in Eq. vanishes 



identically 22]. Therefore, to one-loop order, the contri- 
bution to the EP in the weak field case from the Higgs 
sector is independent of YiJ and is given by Eq. (|5D|) . 

In the weak field limit, the contribution from the 
fermion sector is given by 



J 



Nf 

i—1 n 



(2^ 



E^^E 



2(YiJ)^ 



(27r) 
K + k2 



C^^fc r, r 2 



k^ + m/?] 



(51) 
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where the sum runs over the number of SM fermions, Nf, 
with masses uif^ = -^v- We emphasize that the only 
fermion mass we keep in the analysis is the top mass. 



The first term in Eq. (jSip represents the fermion ideal 
gas contribution [28J , which is explicitly given by ^26| 



2—1 n 



(fik 

(2^ 



Nf 



i=l 



7 

180 



12 



167r2 



■In 



9 

2^2 



Oirrift)}. (52) 



The second term in Eq. (|5ip is subdominant, after taking 
care of renormalization and running of the coupling con- 
stant, as we sketch in the Appendix. Therefore, to one- 
loop order, the contribution to the EP in the weak field 
case from the fermion sector is only given by Eq. ([5^ . 

Finally, since in the symmetric phase the gauge bosons 
do not couple to the external field, their contribution to 
the EP is given by ^ 



(1) 



G n 

IE 



(2^ 



(27r)3 



{61n(^2 



k2 



(53) 



where on the right-hand side in the first line, the index 
G runs over the four SM gauge bosons and thereafter, 
we have used the mass eigenbasis with = g^v^/A 
and m% — {g^ + g'^)w^/4. The factors in front of each 
contribution correspond to the two W's, the Z and the 
photon polarizations. 

Using Eq. ([50)1 into Eq. ([SH)) . the contribution to the 
EP from the gauge boson sector is 



(1) 



gb 



11^^ i2ml- + ml)T^ (2m^. + m|)r ^ mw\ _ ( rnz\ . 4 ^ 



r 



In our expressions for the 1-loop potential, ultraviolet 
divergences are absorbed by the standard renormaliza- 
tion procedure 



X Tr ln[l + n^'99«(o)Z5^(^„,k;m,)], (55) 
for the scalar case, and by 



B. Ring diagrams 

It is well known that the next order correction to the 
EP comes from the so called ring diagrams. These are 
schematically depicted in fig. [5l Their contribution to 
the EP is given by 



^ |E^[-^H(c^„,k;m,)n^'««^(0)] 



N 



4 



i=l 



_fk_ 

(27r)3 



d^k 



TM E]^[-nMA(o)i^,\''(^„,k) 



N 



(56) 



for the gauge boson case. The dominant contribution 
comes from the mode n = in Eqs. ([55]) and dS S]) J28| . 
The scalar case has been treated in detail in Ref. [22] for 
the weak field limit and here we just quote the result 
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V, 



(ring)_ 



Higgs 

4=1 



T 

12^ 



Hi 



V 3/2 



{YHf 

All 



Hi 

48 



T 



(77if + ni)3/2 



(57) 



where the contributions from aU scalars has been ac- 
counted for. 

In order to compute the contribution from the gauge 
bosons, we have to diagonaUze the matrix product 
I¥^^y^{Qi)Dl^ [un, k) in Eq. In the gauge group space, 

we notice that since n°^(0) is, according to Eq. (|44p . di- 
agonal we can use the same matrix that diagonaUzes the 
mass matrix, given exphcitly for instance in Ref. [s^ . On 
the other hand in Lorentz space, the Euchdean version 
of the gauge boson propagator in the Landau gauge can 
be written as 



1 



fc2 



(58) 



h h 

pL ^ c/jy _ ^m"''^ _ pT 



(59) 



It is easy to see that considering only the dominant term 
in the product Uf^\{0)D^'^ {ojn,^) one gets 



ttQ 



1 + 



{k ■ uf 



Ql- (60) 



where 



p'^ 



Considering the n — term and taking the trace, Eq. ((56|) 
gives 



-— {2mUT, H) + TO|(r, H) + nvt^{T, H) - 2™^^ ^ ni%] , 



(61) 



where 'mQ{T,H) are explicitly given by 



{T,H) = 



mliT,H) 



ml{T,H) 



mj, + ig' + g")UG 



[ml + {g^+g'^)UGY 



g^g'^G , 2 {g^ + g'^)r. 1 
\ml + — ^IIg] 



(5' +.9") 



(62) 



Note that both the temperature as well as the magnetic 
field corrections to the gauge boson masses are encoded 
in IIg. When 7? — > the gauge boson masses tend to 
their usual thermal masses 12011 



m|(r,i/^0) ^ m|(r) 
ml{T,H^O) -> m^jT). 



(63) 



where TOg.(T) can be obtained from Eqs. ((44|) and ([62]) 
by setting H = 0. 
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firing) 



To bring about the explicit dependence of y\ 
the hypermagnetic field strength, let us expand the the 
cube of the gauge boson masses appearing in Eq. lpT|) up 



to order (YiJ)^, which is our working order. From the 
expressions in Eqs. (j62p this gives 



(ring) 



(v) 



T 



— {2m3^(T) + m|(r) + (T) - 2 



127r 

T r3 



— {- [2g'mw(T) + [g^ + g'^){mz{T)Rz - m-,{T)R^)\ 




(64) 



where 



Rz = 



2m|(r) 



ml{T) + mUT) 



2ml{T) 



9^9'^ 



mUT) 



{9' 



(65) 



The final expression for the effective potential is ob- 
tained by adding up the results in Eqs. ((46)) . ([50]) . ((52)) . 
dSl, dSZl) and 



V{V) = y,.ee(«) + + yr + n 

, ,(rmg) ,(ririg) 



.(1) 



.(1) 



(66) 



Notice that the dangerous terms which come from 
Lilians cancel with the corresponding terms coming from 



V, 



19) 



Higgs 



The cubic term in the masses of the gauge 
bosons do not cancel but since these masses are posi- 
tive definite, these terms do not pose any problem. Also, 
in order for the terms involving the square of the bosons' 
thermal mass to be real, the temperature must be such 
that 



T>Ti = 



-16ml{v = 0) 
3.g2 -I- -I- 8A + 4/2 



(67) 



which defines a lower bound for the temperature. A more 
restrictive bound is obtained by requiring that yH < fh? 
for the weak field expansion to work. This condition 
translates into the bound 



T > T2 



yH - l&ml{v = 0) 
3.g2 + g'2 + 8A + 4/2 



(68) 



The relevant factor that enhances the order of transi- 
tion, present both in V^^^^] and Vgl'^'^^\ is {yH)'^/rh^ 
which can be traced back to the boson self-energy dia- 
grams involving a tadpole of hypercharged scalars in the 
presence of the external field. 



FIG. 5: Schematic representation of ring diagrams, that con- 
sist in the resummation of successive insertions of self-energies 
in vacuum bubbles. 



VI. SYMMETRY BREAKING 

In order to quantitatively check the effect of the mag- 
netic field during the EWPT, we proceed to plot as 
a function of the vacuum expectation value v. For the 
analysis we use g' = 0.344 and g — 0.637, mz = 91 GeV, 
mw — 80 GeV, / = 1, A = 0.11 which corresponds to 
the current bound on the Higgs mass. 

Figure [S] shows Veff for different temperatures and 
H = 0. This figure shows the usual behavior whereby 
at high temperature the symmetry is restored and, when 
decreasing the temperature, Veff develops a secondary 
minimum that becomes degenerate with the original one 
for a critical temperature T^, where v ^ v'^. Lowering 
further the temperature below produces the system 
to spinodaly decompose. 

To study the effect of the hypermagnetic field on the 
effective potential, we parametrize H — /i (100 GeV)2. 
In what follows, we always choose values for h that com- 
ply with the condition /i(y(100 GeV)2) < ml{v = 0), 
as required from our hierarchy of scales, since the small- 
est possible mass is always that of a scalar boson in the 
symmetric phase. Figure [7] shows Veff for different field 
strengths H3 > H2 > Hi and constant T where T is 
taken as the critical temperature for the h — case. We 
can see that for higher field strengths, the phase tran- 
sition is delayed, favoring higher values of the ratio of 
the vacuum expectation value (v) in the broken symme- 
try phase to the critical temperature Tc- This makes the 
transition stronger first order as the strength of the field 
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0.25 




v/T 



FIG. 6: (Color on line) V^ff as a function of v for differ- 
ent temperatures and H — 0. At high temperature (Ti) 
the symmetry is restored. Decreasing the temperature {T2) 
causes 14// to develop a secondary minimum that becomes 
degenerate with the original one for a critical temperature 
r° < Ti < Ta, where V = V . Lowering further the tempera- 
ture below produces the system to spinodaly decompose. 
For the analysis we use g' = 0.344 and g = 0.637, mz = 91 
GeV, mw = 80 GeV, / = 1, A = 0.11. 



increases. 

Figure [5] shows a comparison of Vef / for different field 
strengths H3 > H2 > Hi at their corresponding crit- 
ical temperatures with TcQit,) < Tc{h2) < Tc{hi) for 
/i3 = 0.25, /i2 = 0.2 and hi = 0.1. From this figure 
we also note that increasing the intensity of the field, 
the barrier between minima becomes higher and the ra- 
tio {v)/T becomes larger, making the EWPT become 
stronger first order @. 

The effect of the hypermagnetic field is therefore 
twofold: first it delays the beginning of the phase transi- 
tion and second the Higgs vacuum expectation value in 
the broken symmetry phase becomes larger. Both effects 
favor the suppression of the sphaleron transition rate and 
therefore help the freezing of a possible baryon asymme- 
try during the EWPT. 



VII. DISCUSSION AND CONCLUSIONS 

In conclusion, we have shown that in the presence of 
an external magnetic field, the EWPT becomes stronger 
first order. Our treatment has been implemented for 
the case of weak fields for the hierarchy of scales viJ ^ 
<C where m is taken as a generic mass involved in 
the calculation. Notice that when this relation is applied 
to the case of the scalar masses, should be regarded as 
m^. We have explicitly worked with the degrees of free- 
dom in the symmetric phase where the external magnetic 
field belongs to the U{1)y group and therefore properly 
receives the name of hypermagnetic field. The calculation 
is carried out up to the contribution of ring diagrams to 
the effective potential at finite temperature. To include 
the effects of this external field, we have made use of the 



0.25 



0.2 - - - h. 




0.1 0.2 0.3 



v/T 

FIG. 7: (Color on line) Veff as a function of v for constant 
T and different values H = /i(100 GeV)^ where hi = 0, 
h2 — 0.03 and /is = 0.06. The choosen temperature T is the 
critical temperature for hi = 0. For higher field strengths, the 
phase transition is delayed, favoring higher values of {v}/Tc. 
For the analysis we use g' — 0.344 and g — 0.637, mz = 91 
GeV, mw = 80 GeV, / = 1, A = 0.11. 



Schwinger proper-time method. In this way, the contri- 
bution from all Landau levels has been accounted for. 
We have carried out a systematic expansion up to order 
(yH)^. The presence of the external hypermagnetic field 
gives rise to terms in the effective potential proportional 
to 1/to^, where m.f — mf -\- Hi, coming from tadpole 
diagrams in the boson self-energies where the loop par- 
ticle is a hypercharged scalar and are their masses. 
These terms are the relevant ones for the strengthening 
of the order of the phase transition since rhi(v) is small 
for small values of v, which enhances the curvature of the 
effective potential near v = 0. 

The presence of the hypermagnetic field has two simul- 
taneous effects: first it delays the beginning of the phase 
transition as compared to the case with no hypermag- 
netic field. Second, the Higgs vacuum expectation value 
in the broken symmetry phase becomes larger. Both ef- 
fects favor the suppression of the sphaleron transition 
rate and therefore help the freezing of a possible baryon 
asymmetry during the EWPT. Although our study sug- 
gests that for reasonable values of the magnetic field 
strength the ratio {v)/T does not reach the necessary 
values to avoid the sphaleron erasure, we should keep in 
mind that we have worked under very restrictive assump- 
tions, derived from our assumed hierarchy of scales. Lift- 
ing some of the restrictions, in particular the weakness 
of the magnetic field compared to the mass scale, could 
open a window to make the above ratio larger, while re- 
maining within the observational bounds. This is a sub- 
ject under current study and will be reported elsewhere. 
The calculation has been carried out in the framework of 
perturbation theory but in order to make contact with 
the values for the masses of the top quark and the cur- 
rent bounds on the Higgs mass, in the numerical analysis 
we have used large values for the coupling constants, in 
particular we have considered the top Yukawa coupling 
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T(h,) 
T(h,) 



H 




v/T 



FIG. 8: (Color on line) V^ff as a function of v for different 
hypermagnetic field strengths /13 > /12 > hi at their cor- 
responding critical temperatures Tdh-i) < Tc{h2) < Tc{h\), 
where H — h (100 GeV)^. We note that increasing the inten- 
sity of the field, the barrier between minima becomes higher 
and that the ratio {v)/T becomes larger. For the analysis we 
use g' = 0.344 and g = 0.637, mz = 91 GeV, mw = 80 GeV, 
/ = 1, A = 0.11. 



f — I. In this sense, our calculation has to be consid- 
ered as an analytical tool to explore this non perturbative 
regime. 

Our results go along the same direction as the ones 
obtained at tree [l3| and one-loop [l3|levels. We should 
however point out that the same problem has been also 
treated in Refs. [18, 19] in the context of strong external 
magnetic/hypermagnetic fields. These authors conclude 
that the magnetic field gives rise to logarithmic terms of 
the ratio of the temperature to the fcrmion masses and 
that for light fermions, these terms increase the inertia of 
the Higgs field producing the phase transition to become 
second order. Similar logarithmic terms appear in a weak 
field expansion of the effective potential in Ref. [131 al- 
though, as pointed out by the authors of that work, their 
weak field expansion is not necessarily reliable since it is 
given in terms of a Borel summable rather than a conver- 
gent series. In this work we have not come across such 
terms but a detailed comparison between the methods 
used in the above mentioned works with ours to include 
the effects of the magnetic field is certainly called for. 
Also the study of the gauge parameter dependence and 
the inclusion of small but not necessarily negligible terms 
proportional to the gauge boson masses is a natural ex- 
tension of this work. Progress in this direction is being 
made and we will report on it elsewhere |35|. 
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Appendix 

First we sketch the computation of the diagram de- 
picted in fig. [3)^//) representing the fermion contribution 
to the one-loop gauge boson self-energy tensor up to or- 
der (yH)^, where it is not evident that the off-diagonal 
components are negligible. We use the notation "II^i,, 
n = 0, 1, 2, where n is the power of yH. To zeroth order 
we have 



(27r)3 
16kuK - Sk'^Sui, 



[k^ + mf'j] [{k — p)2 + nifj 



(69) 



Using the result in Eq. (|29p for fermions we get 



2^2 1 ^f. 
3 2^2 



' ^1 , V 

- + 7_E + In —J 



2 r 



(47r)2 



1 4^2 1 



(70) 



To first order we have 



v«)(f)^E 



d^k 



(27r)3 

X [k^ikf^Ui, + k^Uf^) + ujri{ki_,b,y + k^bi_,) 
- 2cj„/c3(5^j,] [I12 (k, q) + I21 (k, q)] 
+ Sie^uc^p [ksk^u^ + iOnk"b^] 
X [/i2(k,q)-/2i(k,q)]. 

Using again Eq. (P^)) we explicitly get 



(71) 



(yH) 



1 , ifi 

- + 7b + In — 
e TTi 



" V4y 2^2 

X (Sf^sSr^Q+Sf^oSr^s)- (72) 

To second order we get contributions from two terms 



+ 2uJrik3{u^b^ + u^b^) + ujl{2bf_,b^ - 6^^)] 
X /22(k,q), 



(73) 



and 



X. = -8(Yff)^(f)VE / 



d^k 



X {\k^Jz^[kl^uj'i)^2[(k^,b^ + k^b^)k-i 

- (fc^fo^ + fc^6^)c^„)(t^,Hk^)] 
X (/4i(k,q)+/i4(k,q)) 

- 2ie^,,;3 {k^b'^kf' - uj^v^k^\ e 
X (/4i(k,q)-/i4(k,q)). (74) 
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Using once more Eq. (|29)) we get 



{2 + 'Sb'b')S^ 



7^(3) 
487r4 



(75) 



and 



\—) \V 144^ 



(76) 



From Eq. ([TO)) , we see that the leading contribution 
corresponds to the term proportional to and thus 
the result in the first of Eqs. (|42p. after summing over 
all hypercharged fermions. After factorizing the leading 
contribution, we also notice that the terms proportional 
to (mf^/T)'^ come together with a factor reminiscent of 
renormalization and running of the fermion mass and 
couplings with the scale fi. We do not carry out this pro- 
cedure in detail but when absorbing these factors into the 
redefinition of the fermion mass, we see that corrections 
to the leading term are proportional to (m^^/T)^ and 
therefore for the values taken by v during the EWPT, 
these corrections can be safely neglected. 

From Eq. l[7^ . and by the same procedure as for the 
case of the zeroth order terms, the off-diagonal terms 



contribute an amount proportional to Y_ff/T^, compared 
to the leading term and thus for the hierarchy of scales 
considered, this term can also be neglected. 

Finally from Eqs. ([75)1 and ([75]) we can see that the 
second order contributions are proportional to {yH/T'^Y 
and therefore can also be safely neglected. 

Next, we sketch the computation of the hypermagnetic 
field contribution to the fermion one- loop effective poten- 
tial V^'^\ According to Eq. ([5T|) . this contribution is of 
second order in yH. 



(1) 



(2^)3 (cj2 +k2+m/2)3 



(77) 



Using Eq. ^ we get 



'■V. 



(1) 



2(YiJ)2 ri 7^(3) 



327r2 



27r2 r2 



ln(^)+7.+ 3 



(78) 



After renormalization and running of the couplings with 
the scale we see that this contribution is of order 
{yH/T^Y ^-nd can also be safely neglected. 
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